Collocation methods for the solution of von-K arm an dynamic non-linear plate systems

The von-K arm an nonlinear, dynamic, partial differential system over rectangular domains is considered, and numerically solved using both the Chebyshev-collocation and Legendre-collocation methods for the spatial discretization and the implicit Newmark-b scheme combined with a non-linear fixed point algorithm for the temporal discretization. As the system is non-linear, involving operators of different orders, different timescales, and may contain initial/boundary incompatible conditions at the domain’s corners, it is our aim to highlight some of the difficulties inherent in its numerical treatment using collocation methods. We begin by examining the Chebyshevcollocation scheme considered in [Comput. Meth. Appl. Mech. Engrg. 193 (6–8) (2004) 575]. In that work, filtering was used to stabilize the full von-K arm an system, however the source of the instability was not discussed. We first show that IC/BC incompatibilities are not the culprits and demonstrate empirically that the source of the instability was the numerical treatment of the linear cross-derivative terms found in the in-plane equations. We prove mathematically that although the continuous linear system is well-posed, the Chebyshev-collocation solution of problems involving cross-derivatives and second-order time derivatives is unstable. Instead of adding sufficient filtering to counter the energy growth due to the linear terms, we employ the Legendre-collocation method. We show empirically and prove that the Legendre-collocation scheme for the linear system has the necessary characteristics to remain stable. Nevertheless, the collocation treatment of the non-linear terms still causes the solution to be unstable at long time. We thus return to advocating a filtering solution, but in this case only filtering the nonlinear terms to compensate for the inherent aliasing error due to our collocation treatment of these terms. We conclude by comparing the results of filtered Chebyshev-collocation solution, non-linear terms filtered LegendreCorresponding author. E-mail addresses: [email protected] (Z. Yosibash), [email protected] (R.M. Kirby), [email protected] (D. Gottlieb). 1 Research performed while the author was on a Sabbatical leave at the Division of Applied Mathematics, Brown University, Providence, RI, USA. 0021-9991/$ see front matter 2004 Elsevier Inc. All rights reserved. doi:10.1016/j.jcp.2004.03.018 Z. Yosibash et al. / Journal of Computational Physics 200 (2004) 432–461 433 collocation solution, and a simplified von-K arm an (solution omitting the inertial terms in the in-plane equations) as presented in (loc. cit). 2004 Elsevier Inc. All rights reserved.